Decoherence of Einstein-Podolsky-Rosen pairs in a noisy Andreev entangler 
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We investigate quantum noise effect on the transportation of nonlocal Cooper pairs accross the 
realistic Andreev entangler which consists of an s-wave superconductor coupled to two small quantum 
dots at resonance which themselves are coupled to normal leads. The noise emerges due to voltage 
fluctuations felt by the electrons residing on the two dots as a result of the finite resistances in the 
gate leads or of any resistive lead capacitively coupled to the dots. In the ideal noiseless case, the 
setup provides a trustable source of mobile and nonlocal spin-entangled electrons and the transport is 
dominated by a two-particle Breit-Wigner resonance that allows the injection of two spin-entangled 
electrons into different leads at the same energy [P. Recher, E. V. Sukhorukov, and D. Loss, Phys. 
Rev. B 63, 165314 (2001)]. We seek to revisit the transport of those nonlocal Cooper pairs as well 
as the efficiency of such an Andreev entangler when including the quantum noise (decoherence). 

PACS numbers: 73.23.Hk, 73.40.Gk, 73.50.Td 
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The concept of nonlocal pairwise-entangled quan- 
tum statesA, namely the Einstein-Podolsky-Rosen (EPR) 
pairs^ is extremely fundamental for testing the violation 
of Bell inequality 3 and at the same time is also cru- 
cial for efficient quantum communication* and quantum 
teleportatior£ for example. Tests on the entanglement 
of massless particles like the photons already exists but 
not yet for massive particles such as the electrons. A con- 
crete challenge is to realize an entangler of electrons, i.e., 
a device that generates spin singlets that are made out 
of two electrons. Recently, an interesting setup involving 
an s-wave superconductor weakly-coupled to two quan- 
tum dots, which themselves are coupled to normal leads 
has been envisioned in Ref. 0- The spin correlations of 
an s-wave superconductor induce a spin-singlet state be- 
tween two electrons, each of which can now reside on a 
separate quantum dot; this results in the formation of a 
nonlocal Cooper pair described by the quantum state 



\dd) = [dt ?< 4 - 44] 



(1) 



where dl produces an electron with spin a on the dot I 
and \i) = |0)s|0)d|^) (|Q)s is the quasiparticle vacuum 
for the superconductor, |0)rj means that both dot levels 
ei are unoccupied, and defines the occupation of the 
leads which are filled with electrons up to the electro- 
chemical potential fj,i) is the initial state. A prerequisite 
to perform such a nonlocal spin-entangled electron state 
is the Coulomb blockade phenomenon which typically for- 
bids double occupancy on each quantum dot: the spin 
singlet coming from the Cooper pair remains preserved 
in this process even though the two involved partners 
are well separated physically because they reside on dif- 
ferent quantum dots. More precisely, such an Andreev 
entangler— of EPR pairs has been predicted to be well ef- 
ficient wheni £ ft > (k F 5r) where = 1/(ttA) + 1/U, 
U being twice the charging energy of a given dot, A the 
superconducting gap, Sr = \Sr\ which might be of the 
order of the distance between dots denotes the distance 
between the points on the superconductor from which 
electrons 1 and 2 tunnel into the dots, hp (kp — \kp\) 



the Fermi momentum in the superconductor, and 7 the 
effective tunneling rate between the dots and the normal 
leads. We have assumed a two-dimensional supercon- 
ductor. The parasitic direct Andreev process where two 
spin-entangled electrons from the superconductor navi- 
gate accross the same dot can be minimized due to the 
Coulomb blockade effect. In principle, assuming that the 
tunnel barriers between the dots and the superconduc- 
tor are weak enough, this allows to suppress the other 
parasitic quantum process namely the superconducting 
cotunneling^ (transfer of an electron from one dot to an- 
other via the superconductor). In this setup, there is a 
direct correspondance between the idea of preserving the 
spin entanglement on the two dots - implying that the 
distance between the two dots is smaller than the coher- 
ence length £ of the superconductor - and the result to 
get a two-particle Breit-Wigner resonance allowing the 
spin injection of two spin-entangled electrons into differ- 
ent leads at exactly the same orbital energy and hence 
a quite prominent stationary current when a bias volt- 
age is applied between the superconductor and the two 
leads. Remember the one to one correspondance between 
non-local spin-entanglcmcnt on the two dots and promi- 
nent stationary current between the superconductor and 
the leads triggered by the crossed Andreev reflection. A 
similar conclusion arises when replacing the single-level 
quantum dots by one-dimensional Luttingcr liquidsASiii 
however the efficiency condition of the entangler is mod- 
ified asi^ (A/2/i) 27 '' 3> (kpSr) where /i is the difference 
of electrochemical potentials between the superconductor 
and the Luttinger leads and j p = (K p +K~ 1 )/4 : -l/2 > 0, 
where K p < 1 is the Luttinger parameter, stands for the 
exponent for tunneling into the bulk of Luttinger leads. 

The nonlocal spin entanglement might be detected 
through (current) noise measurements^ 3 -. In fact, using a 
beam splitter in the setup of Ref. Q, Samuelsson, Sukho- 
rukov, and Biittiker have shown that the electrical cur- 
rent noise could serve as a test of spin-entanglement^. 
Following Ref. Ha. a diagnosis on spin current noise and 
Bell inequalities in the same setting has been done by 
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FIG. 1: (Color online) The noisy entangler: Two spin- 
entangled electrons forming a Cooper pair tunnel with ampli- 
tude Tsd from points ri and r-2 of the superconductor (SC) to 
two dots D\ and D2 by means of crossed Andreev reflection. 
The dots are coupled to normal leads L\ and L2 with tunnel- 
ing amplitude \Tdl\ > |Ts£>|- The SC and leads are main- 
tained at chemical potentials (is and fii respectively so that 
/1 = /is — m > 0. We show the dissipative leads (which are 
not necessarily the gate leads) embodied by the impedances 
Z\ and Zi producing voltage fluctuations (8Vi(t)) on each dot. 



Sauret, Martin, and Feinbergi&. Bell-inequality checks 
in solid-state systems can be thought as interesting gen- 
eralizations of the corresponding tests with photonsi£*i&. 

Many theoretical papers are devoted to the crossed 
Andreev reflection (the elastic emission of spin-entangled 
electrons in different leads) as well as to the distance- 
contribution of the crossed Andreev reflection^ and to 
its implications on the conductance matri:s22i2ii22i22i or 
on the current-current cross correlation^. 

On the other hand, the destruction of quantum- 
mechanical phase coherence due to coupling of a sys- 
tem to a dissipative bath is a subject of great impor- 
tance partly because of its connection to fundamental 
issues related to the quantum measurement process or 
the quantum-classical crossover. Some specific condensed 
matter setups have been also proposed to answer the 
question whether decoherence at zero temperature is pos- 
sible at allSSiS 6 ^. Following this line of ideas, it is then 
interesting to ask to what extent zero-point fluctuations 
of the electrical environments of the two dots will affect 
the non-local Cooper pair of charge 2e and hence the sta- 
tionary current(s) discussed in Ref. 0- Therefore, we now 
turn our attention to the specific setup of Fig. I where 
the external impedances Z\(w) and Zi(u>) show resistive 
behavior at low frequencies. Albeit voltage fluctuations 
from gate leads in GaAs, in double dot geometries, are 
thought to be quite smali2& and often ignored, the quan- 



tum noise may emerge from any resistive source capac- 
itively coupled to the dots. More precisely, we have in 
mind the GaAs heterostructures of Ref. [23 where quan- 
tum dots are built in the proximity of two-dimensional 
electron gases with (tunable) resistances which can even 
reach few kQ when exploring the low density limit of 
those two-dimensional electron gases. The noise might 
also stem from the metallic leads LI and L2 on Fig. 1; 
indeed, similar results can be derived. The quantum 
noise will couple to the nonlocal Cooper pair through 
the requirement that to hop on the dot an electron must 
transmit a part of its energy to the electrical environment 
(bath) and thus the \DD) state must be thought as: 

\DD)\D B ) = e-^idl^ - dl^MiB), (2) 

describing explicitly the "entanglement" of the nonlo- 
cal Cooper pair to the quantum noise; \%b) embodies 
the initial bath (ground) state, \D B ) = e~ iS ^\iB) de- 
scribes the bath excited state when the two electrons 
reside on their respective dots and at a general level 
5(j){t) = {e/h) J^(SVi(t r ) + SV 2 (t'))dt' embodies the phase 
conjugate to the sum of voltage fluctuations felt by 
the two dots^i. Note that Eq. (2) stems from the 
fact that in the presence of voltage fluctuations, the 
dot's electron annihilation operators will be modified 

as2£ di a — > di a e^ e ^ n ^o SVl ^ ' dt . As a matter of fact, 
in the case of moderate ohmic resistors i?i and R2 such 
as Zi{lo) = R t , (e-^W) = {i B \e~ iS ^\i B ) always goes 
to zero at sufficiently long times both in the finite tem- 
perature regime and in the quantum realm2& (at finite 
temperatures, one gets a relatively fast exponential de- 
cay versus time whereas approaching the quantum limit 
one reaches a less severe "power-law" decay) . By averag- 
ing over the bath degrees of freedom, we can rewrite Eq. 

(2) as \DD)(i B \D B ) = (e"^) [c4 T 4i ~ hence 
we infer that the environments will contribute to the de- 
struction of the (nonlocal) Cooper pair object in the long- 
time limit. Our interest is to study the consequence of 
the nonlocal Cooper pair long-time "decoherence" on the 
current between the superconductor and the leads. By 
decoherence we explicitly refer to the long-time decay of 
the noisy Cooper pair through the factor (e^ 5 ^*)) even 
though the spin entanglement is not destroyed. 

Even though electron spins in semiconductor het- 
erostructures are in principle very robust to decoherence 
with dephasing time approaching microseconds^, in the 
realistic setup of Ref. one can nevertheless discuss 
the phenomenon of decoherence (decay) of a nonlocal 
charged-2e Cooper pair during the transportation pro- 
cess due to the electrical circuits existing in the vicinity 
of the dots; this is the subject which will be addressed in 
this paper. Our paper is organized as follows. In Sec. I, 
we introduce the model of the noisy Andreev entangler 
and subsequently we write a general formula for the cur- 
rent resorting to the T-matrix approach. In Sec. II, we 
discuss the result for a single environment — assuming 
for example that R2 = — both in the limit of large resis- 



tances R\ /Rk ^> 1, Rk = h/e 2 = 25. 8kfl being typically 
the quantum of resistance, and in the maybe more real- 
istic situation of moderate resistances R\/Rk « 1. In 
Sec. Ill, we extend our results for the two bath situation 
where both R\ and i?2 are assumed to be non-negligible. 
In Sec. IV, we also take into account the parasitic direct 
Andreev processes where the two electrons issued from 
the superconductor eventually jump onto the same dot; 
this processes are nevertheless reduced at low temper- 
atures due to the prolific combination of the Coulomb 
blockade principle which typically forbids two electrons 
to jump simultaneously onto the same dot and the su- 
perconducting gap which avoids quasiparticle formation 
in the superconductor. In the conclusion, we seek to re- 
visit the condition of efficiency for such a noisy entangler. 
Appendices are devoted to the details of calculations. 



I. MODEL OF THE NOISY ENTANGLER 

A. Hamiltonian of the noisy Andreev entangler 

Similar to Ref. 0, we exploit a tunneling Hamiltonian 
description of the system, H = Hq + Ht + H^ th , where 



^ y jjbath 
I 



(3) 



The superconductor is embodied by the usual BCS 
theory^ which for convenience has been summarized 
in Appendix A. In brief, H s = Y,ka ^k7k<j7ko- where 
a =T,J, represents the spin index, 7k CT describe exci- 
tations out of the BCS ground state |0)s defined by 
7ko-|0)s = 0, i?k = vA^+^k i s * ne quasiparticle en- 
ergy, and £k = £k — Ms is the normal state single-electron 
energy counted from the level us where live the Cooper 
spin-singlet particles. Both dots are embodied by a single 
level with energy t\ very close to /is and are typically gov- 
erned by an Anderson model Hd 1 = Q Y) c d, dicr+Un^ni 
and I = 1,2. The resonant dot level e; can be adjusted by 
the related gate voltage (or by Vi on Fig. 1). Other levels 
do not participate in the transport when the level spacing 
of the dots is sufficiently large implying Se > fi > ksT. 
Again, fi = /is — fJ-i is the difference of electrochemi- 
cal potentials between the superconductor and the leads. 
Moreover, through the on-site Coulomb U repulsion a 
double occupied state is rather hindered to form on each 
dot; U is equal to 2E C where E c = e 2 /(2C) is typically 
the charging energy on each dot and C denotes the to- 
tal dot's capacitance. Keep in mind that this Coulomb 
blockade argument stands for a key point in the efficiency 
of this Andreev entangler of EPR pairsi The leads are 
normal and embodied by a non-interacting theory Ht, = 
Skcr £ ka/kcr a 'kcr- We have to consider final two-particle 

states of the form |/) = (l/>/2) [a{ pT 4 qi - 4pi a 2qT]l*) 
with energy ej = e p + e q . The preserved spin singlet 
state is formed out of two electrons, one being in the p 
state in lead 1 while the other one is in the q state in 
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FIG. 2: (Color online) Impedances as transmission lines. 



lead 2. Since the total spin is conserved, the singlet state 
of the initial Cooper pair will be conserved in the trans- 
port process and the final state must satisfy S z = 0. The 
S z = configuration of the triplet state is excluded as 
long as the distance between the two dots is smaller than 
£. Tunneling events from the dot I to the lead / or to the 
point r; in the s-wave superconductor is described by the 
tunnel Hamiltonian Ht — H$d + (Hdl + Hdl)> with 

H S d = T SDdlMn) + h.c. = H SD t , (4) 

la I 

Hdl = ^2T DL a f lka di a = ^2H DlLl . 

Ika I 

Here, ip a {ri) annihilates an electron in the superconduc- 
tor at the site r; and d\ creates it again on dot I with 
amplitude Tsd- Tunneling from the dot (lead) I to the 
lead (dot) I is described by the tunneling amplitude Tjy^ 
{T DL ). The k dependence of Tdl can be safely neglected. 
Moreover, like in Ref. 3, we require the dot-lead cou- 
pling to be much stronger than the superconductor-dot 
coupling, i.e., \Tsd\ < \Tdl\, so that electrons that enter 
the dots from the superconductor will leave the quantum 
dots to the leads much faster that new electrons can be 
provided to the dots from the superconductor. Addition- 
ally, a stationary occupation due to the couplings to the 
leads is exponentially small if fi > ksT. Thus, in the 
asymmetric barrier case, the resonant dot levels e/ are 
occupied only during a virtual process. Other important 
parameters are the tunneling rates 7; = 2hvi\Tdl\ 2 and 
75 = 2ttvs\Tds\ 2 where v\ is the density of states per spin 
of the leads at the Fermi energy and vg for the supercon- 
ductor will be defined as 1/A. Recall that we will work 
in the regime where 7; > 75 and A, U, Se ^ fj, > 7/, UbT 
and close to the resonant condition for the dots ei « [is- 
We model the impedance Zi(ui) in a microscopic fash- 
ion through a long dissipative transmission line composed 
of an infinite collection of L t iCu oscillators (Fig. 2). Our 
environments are modeled in a usual way akin to Refs. [2^ 
an d33: the charge operator Q n i on the capacitor between 
two inductances L t i and the conjugate flux (Ti/e)<j) n i are 
mapped onto the operators^ Qi(x) and 4>i{x) which are 
precisely described by the diagonalized Hamiltonian 



H, 



bath 



2C tl 



-— sin^ — Uf(x) . (5) 



The charge (fluctuation) operator Qi(x) and the 
phase operator 4>i(x) obey the commutation relation 



4 



[4>i (x) ,Qi(y)/e] = iS(x — y). The Hamiltonian containing 
the couplings with the dots reads 



rrbath 

H D = 



la 



-prQoi4 a di a = y^eSVid^dia. 



la 



(6) 



This term may arise from the extra capacitive cou- 
pling between each dot and the voltage fluctuations (the 
quantum noise) SVi(t) = Qoi/Ci with Qqi denoting the 
charge fluctuation operator on the given capacitor C/, 
emerging from the finite impedance Zi(uj)-^2&. Ac- 
cording to Ref. one can thoroughly identify Qqi = 
\[2 J dxcos(nx/2)Qi(x). At low frequency u> < uj c i = 

\/{RiCti) where Ri = \J Lti/Cti, the transmission line 
provides an impedance Zi(oj) = Ri/(1 + itu/uj c i) « R[. 

Below, we will absorb the H h £ th coupling into the tun- 
neling terms through the unitary transformatio: 



.31 



U = exp 



iJ2Wi(t)dld, 



(J) 



where we have defined 8<j>i{t) = (e/h) f dt'5Vi(t') = 
v / 2/ o L dxcos(7rx/2)0 / (x). We get Ud la W = e i5 ^ d la 



and Ud\ a W 



e lS( Pi( t )dj a . Moreover, exploiting the 
correspondancc H' = UHU^ + ifr^U^ wc realize that 
the couplings of the dots to the electrical baths can be 
completely absorbed in a redefinition of the tunneling 



Hamiltonian as H' T 



h 'sd 



h 'dl 



H'£ L where 



H'sd =J2 T s° d le- iSMt) Mn) + h.c. = ^H' SDl (8) 



la 



h 'dl - E T DLa\^ a dia<. 



iS<t> t (t) 



Ika l 



At a very general level, the effect of the environments 
can be embodied by a fluctuating phase bound to the 
dot's electron creation and annihilation operators such 
that the total Hamiltonian turns into H' — Hq + H' T . 



B. T-matrix and general current formula 

In the quantum (zero-temperature) regime the current 
of two electrons passing from the superconductor via vir- 
tual dot states to the leads is formally given by 



fi 



(9) 



P.Q 



2c 



Y,27rp i \{f B \{f\T(e i + E B )\i)\i B ) 



x S (e f ~ e> ~ E B + E B ) 
2e x - 

= -T }_^ 27r Pi W i,DDW DD j-, 



where W/i embodies the transition rate from the super- 
conductor to the leads taking into account transitions 



into the electrical environments; (f B \ denotes the final 
(excited) state of the baths when the injected electrons 
arrive in the leads, (E B — E B ) represents the energy sup- 
plied to the environments during the EPR transportation 
process, Wddj and Wi^dd stand for the transition rates 
from the dots to the leads and from the superconductor to 
the dots respectively in the presence of the fluctuations 
in the gate voltages, and pi is the stationary occupa- 
tion probability for the entire system to be in the initial 
ground state where as introduced in the introduc- 

tion \i) — \0)s\0)d\m) and \i B ) depicts the initial state 
for the environments. Along the lines of Ref. 0, to cal- 
culate the transition rates Wddj and Wi^dd we resort 
to the T-matrix approach. The on-shell transmission or 
T-matrix at the energy E is precisely defined as 



T(E) = H' T G{E) 



G (E) 



H'rp 



1 



E + irj- H 



-,{E-H Q ); 



( 10 ) 

we have introduced the Lipmann-Schwinger operators 
G{E) = l/(E- H 1 ) and G = l/(E - H ) as well as 
the small positive real number 77 that we take to zero at 
the end of the calculation. We rewrite the T-matrix as 



T(E) = H' T 



OO 

E 

n=0 



1 



E — H< 



+ if] 



(11) 



It is appropriate to decompose T(E — + E B ) into the 
partial T matrices T' and T" . When two spin-entangled 
electrons from the superconductor leave to the dots 



T" = 



1 



-m 



1 



irj + E-Ho'^^'iv + E-Hn 
+ (l«->2) 



-m 



sd 2 



(12) 



l11 + E -H SDl i V + E-H USD2e 

+ (1^2), 

where (1 «-> 2) refers to the same term but exchanging the 
roles of the labels 1 and 2. T" refers to the (dissipative) 
crossed Andreev process. For the resonant dot <-> lead 
tunneling we must keep all the terms of the series 



= H' DlLl 



(13) 



E 



TTl* TTl 

E-Ho + ir] DlLl E-H + ir] DlLl 



1 



E-H + iri 



■H'n 



E 



1 



Tjl* 

- H DL- 



1 



u] + E-H ^ir] + E-H 



-H' T 



DL 



°° / 1 1 

H 'd 2 l 2 E {e-H q + k, H '^ L * E-H + 1^'°^ 



H'n 



E-Ho + ir] DlL 



E 

m=0 



1 



TTl* 

- H DL- 



1 



E-irj-Ho E — ir/ — Hq 



-H' T 



DL 



■5 




FIG. 3: (Color online) Energy conservation in the absence of 
noise: 2fis — ei + £2 = e P + e q - When an electron on the dot 
1 tunnels onto the lead 1(1) this induces a shift of the Fermi 
energy — > e c (2); e c depicts the edge of the conduction 
band. In the text: e± ~ £2 ~ Ms, Mi ~ A*2 5 an d e c ~ ei ~ £2- 

To compute the resonant dot <-> lead tunneling, we 
have explicitly taken into account virtual \DD) — > l-D-D) 
transitions via the sequences |.DD) — * — > ... — > 
|DD) or |DD) -> \Dq) -> ... -> |DD). Again the 
baths are insensitive to those virtual transitions which 
take place in a very short time and let the state of 
the electron system unchanged. Here, \pD) stands for 
a\ V(J d\_ a \i) and implies that the electron 1 is in lead 1 
(LI) whereas the other electron resides on dot 2. We have 

also introduced the state \Dq) = d\a a \< i -a\i') represent- 
ing the state with one electron on dot 1 and the other 
one in lead 2 (L2). At this point we emphasize that vir- 
tual states with both electrons in the leads leading to 
\DD) — > \pD) — > \pq) — > \Dq) — ► \DD) are suppressed by 
a factor 7;//i < 1 compared to that with only one elec- 
tron in the leads and therefore can be neglected^. Recall 
that \pq) stands for aip<r a 2q-o-K)i where p is the momen- 
tum from lead 1 and q from lead 2; this is the final state 
I/) of the electrons on the leads with energy et = £ q + £ p . 



C. "1-photon" approximation 

In Eq. (13), we have assumed that the tunneling 
events are almost instantaneous implying that only one 
"photon" is emitted in each bath during the two-particle 
Breit-Wigner resonance between dots and leads. In all 
the products Hp L (t)H' D L (0) appearing in Eq. (13) we 
have identified ( e -«&(*) e «<fc(o)) a ^^Ri/Rk ^ ^ 
that is ensured when the time t for an electron on dot I 
to virtually jump in lead I and then to go back to dot I is 
much shorter than l/u> c i ~ RiC t i- Assuming that C t i is 
large enough this should be well satisfied even for weak 
resistances Ri . Note that configurations with emission of 
multiple "photons" in the same bath would only result in 
small corrections in the final current. Keeping only the 
"1-photon" contribution, T may be summarized as^ 



where T' yields the same form as T if one replaces H' D L 
by i.e., without dissipation — R2 — 0); this 

part is at the origin of the two-particle Breit-Wigner res- 
onance between the dots and the leadsi Since \Tsd\ < 
\T DL \ we may always rewrite T" = e -»(^i+^)y". 

Below, we will distinguish between the case of a single 
dissipative bath implying R% = and therefore 8(f>2 = 
and that of two (independent) baths. When R2 = 0, the 
transition rates Wddj and Wi^dd are given by 

Wddj = \{ P q\T'{e i )\DD){f B \e u ^\D B )\ 2 (15) 

x 6(e f - e x - e 2 + E f B - E%), 

W hDD = \{DD\T"(e i )\i){D B \e- iS ' t ' l \i B }\ 2 

x 6(ei + £2 + E B — E* B — e,). 

Note that we have replaced (pq\T'(E)\DD) and 
(DD\T"(E)\i) by (pq\T'(ei)\DD) and (DD\T"(ei)\i) re- 
spectively; this will be thoroughly justified in Sec. II B. 
Furthermore, we have introduced the energy of the inter- 
mediate state of the bath Eg . The energy-conserving S 
functions traduce the fact that for each tunneling process 
the energy of the full system including the bath is con- 
served or that the energy supplied to the bath is equal 
to the energy lost by the electrons. The product of those 
two 5 functions is equivalent to the S function in Eq. (9) 
(dimension of current will be implicitly respected below) . 
When i?i = i?2 = 0, we recover the formulas of Ref. 

W DDJ = \(pq\T'\DD)\ 2 d(e f -e 1 -e 2 ), (16) 
W hDD = \{DD\T"\i)\ 2 5{e x + e 2 -ei), 

or more precisely we obtain the following expression^, 
W fi = 2n \{pq\T'\DD)\ 2 \{DD\T"\i)f S(e, - e t ). (17) 

The chemical potentials t\ and 62 of the quantum dots 
can be tuned by external gate voltages such that the 
coherent tunneling of two electrons into different leads 
is at resonance, described by a two-particle Breit-Wigner 
resonance peaked at £1 + 62 = 2//g = e q + e p (Fig. [3J|; the 
situation we consider is when the two dots are close to the 
resonance condition ei = £2 = 0. We will choose energies 
such that £j = 2/is = 0. It is then an interesting question 
to understand how the quantum dissipation affects the 
elastic tunneling of those EPR pairs. In the case of two 
independent baths, we straightforwardly generalize 

W hDD = \{DD\T"\i)\ 2 \{D Bl \e- i5 ^\i Bl }\ 2 x (18) 

I (D B2 \e~^ \i B2 )\ 2 S(e 1 + e 2 + £ - £ E Bl - e t ), 

1 1 

and similarly 

Wddj = \( P q\T'\DD)\ 2 \(f Bl \e i5 ^\D Bl )\ 2 x (19) 

I {f B2 \e«+* \D B2 )\ 2 x S(e f - e x - e 2 + J2 E k ~ E^)- 

/ 1 



We have decomposed \Db) = \Db 1 )\Db 2 },--- . 

Before to pursue, we shall discuss what is the value 
of pi from Eq. Q in the presence of the environments. 
Similar to the noiseless case, we estimate pi = 1 — 0(7) — > 
1 where 7 = 71 + 72 . More precisely, the initial (gound) 
state \i) — |0)s|0)d|/u/) is such that the highest level 
of the dots is unoccupied, there is no quasiparticle on 
the superconductor which is immediately fulfilled when 
kgT -C A, and the Fermi level of the leads remain fixed 
to p\. Since we consider asymmetric barriers |Tdl| > 
the most prominent transfer of electrons is between 
the leads and the dots. Therefore, the probability for 
the system to remain in the state |i)|is 1 )|is 2 ) after a 
time t can be estimated as 1 — 0(7)- We argue that the 
environments will not affect this equality because as long 
there is no electron on the dots the environments will 
unambiguously remain in their ground states \iBi)- 



II. DISSIPATION ON DOT 1 
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FIG. 4: (Color online) Function Pi(E) for different values of 
a — ocx — Rk / R\ from to 00. When Ri is small one observe 
a blatant singularity at E = converging to a 5(E) function 
when Ri = whereas for huge Ri (ai — > 0) the capacitance 
Ci plays a crucial role leading to Pi(E) — S(E — E c i). We 
analyze the two limiting cases of small and very large Ri. 



Let us assume that R2 — 0. To calculate the rate for 
electron tunneling from the superconductor to the dots 
we have to evaluate W^dd from Eq. (15). Similar to 
Refs. I2I31I we trace out environmental states leading to 



/+00 
■ 
-OO 



dt 



,i{e%— ex— e 2 )t/H 



(20) 



x /gti^i(t) e -i6<j>i(0)\ 



The brackets denote an average over the initial bath 
ground state For later convenience we like to intro- 
duce the abbreviation J x (t) = ([5<f>i(t) ~ (50i(O)]<5^>i(O)) 
as well as the Fourier transform 



dt exp 



Ji(t) + -Et 



This already permits us to write down 



W i<DD = \(DD\T"\i)\ 2 P 1 (e l -e 1 -e 2 ). 



In a similar way, we extrac 
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W DDJ = \(pq\T'\DD)f Pi(ei + e 2 - e f ). 



(21) 



(22) 



(23) 



A. General discussion on P 1 (E) 

We may interpret P\(E) as the probability to emit 
the energy E to the electrical circuit when transferring 
an electron from the superconductor to the dot 1 or from 
the latter to the corresponding lead 1. It is certainly 
useful to know more about the function Pi(E) as well 
as J\(t). We will first envision the case of a very large 
resistance R\. In that case the dissipative lead yields 
an effective impedance Z e ff(ui) = l/(i?^ 1 +iu>C L ) which 
tends to (ir /Ci)5(lu), Ci being the capacitance between 



the dissipative lead and the dot 1. For the correlation 
function J\ (t) this concentration of environmental modes 
at low frequency means that the short-time expansion 31 



J x {t) = 



+ OC 



dw $teZ eff (u;) c _ iult 
uj Rk 



C x Rk 



works for all times. This results in 

Pi(E) =S(E-E cl ), 



it, (24) 



(25) 



so that in order to hop onto the dot 1 an electron must 
transfer to the environment an amount of energy corre- 
sponding to the charging energy E cl — e 2 /(2Ci) of the 
capacitor C\\ this will fatally lead to a Coulomb gap 31 
in the current for low applied bias voltage p <C E c \ be- 
tween the superconductor and the leads (LI and L2). For 
small resistances a\ = Rk/Ri ^ 1 the gate lead may 
be described by the frequency-independent impedance 
Z e ff = Z\ = Ri. Based on the transmission line rep- 
resentation for the environment then we may identify 37 



P 1 {E) 



exp(-27 e /ai) 1 
r(2/ai) E 



7T E 

oti Ed 



2/ai 



(26) 



where 7 e = 0.577... is the Euler constant. The factor 
appearing may be motivated by the behavior of the cor- 
relation function J\ (t) for large times 3 ^ 



Ji(i) 



2 r 

Oil 



ln(aiE c it/irh) + «- + 7e 



(27) 



The function P\ (E) has been summarized through Fig. 
4 with the two distinct behaviors at low and large R\ . 



B. Discussion on tunneling matrix elements 

Now we want to properly justify the fact that the baths 
"cancel out" in the computation of the tunneling matrix 
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elements in Eq. (15). In the limit of a weak resistance 
Ri this is straightforward since the bath 1 only absorbs 
a small amount of energy during the tunneling events 
(Pi(.E) is strongly diverging at E = 0) and therefore in 
all the Lipmann-Schwinger operators appearing in Eqs. 
(12) and (13), for a given bath state \olb), one can al- 
ways formally replace {olb\ J2i Hi ath \otB) — > E l B . Inter- 
estingly, we like to emphasize that for large resistances 
this argument still holds. In the large resistance limit, 
one must typically satisfy Eg — > E B + E c \. Thus, in or- 
der to get a finite current, one must thoroughly re-adjust 
the chemical potential of the SC lead such as £j — > ei+E c i 
with Ci = 0; see Eq. (25). When focussing on the tunnel- 
ing of a Cooper pair from the SC to the dots, hence one 
requires to evaluate (D B \ (DD\ E i +Ec *_ Ho+in ■ Applying 

(D B \H = E° we check that (D B \(DD\ E . g+E ^_ Ho+iri 
is equivalent to (DD\ _ H ] +i by setting H\ ath = in 
Hq. We can thus substitute Go by its expression in 
the absence of the bath 1. This procedure can be ex- 
tended to two baths. For sake of clarity, calculations of 
(pq\T'(ei)\DD) and (DD\T"(e l )\i) (with the substitution 
jjbath _ q j n are derived in Appendices B and C. 

C. Large resistance limit 




FIG. 5: (Color online) At large Ri, the capacitance Ci must 
be taken into account leading to an additional Coulomb gap ~ 
E c \\ to compensate for this gap, one must thoroughly adjust 
Us — > /is + Ed/2 and — > fxi — E c \/2. We will assume 
E c i < A to prevent any quasiparticle poisoning from the SC. 

suppressed by a factor (^/E c i) 2 . This stems from the 
physical fact that shifting the electrochemical potentials 
of the leads such that /i/ — > fii — E c \/2 hampers the 
two-particle Breit-Wigner resonance between the dots at 
resonance (ei + 62 = 0) and the leads. In brief, the appli- 
cation of a prominent bias voltage between the dots and 
the leads somehow produces the "decoherence" of the 
EPR pair hence affecting the crossed Andreev current. 



If we maintain the electrochemical potentials of the 
leads LI and L2 and of the superconductor so that t\ + 
e 2 = 2/is = e p + e q = due to Eq. I|25|) we immediately 
infer that the current will inevitably go to zero. Indeed 

W hDD = \(DD\T"\i)\ 2 5(2^ s -e 1 -e 2 -E cl ) , (28) 
W DDJ = \(pq\T'\DD)\ 2 S(e! + e 2 - e f - E cl ) , 

reflecting the dynamical Coulomb blockade phenomenon 
resulting from large impedances^. On the other hand, 
one could envision to symmetrically modify the elec- 
trochemical potentials of the leads LI and L2 so that 
Ms -* Ms + Ea/2 and fii — ► - E cl /2 (Fig. 5). In that 
case, the two 5 functions above would be satisfied. Tun- 
neling becomes possible only if the energy at disposal is 
equal to E c i. We assume for this circumstance that the 
superconducting gap is large enough A E c \ so that the 
superconductor is not subject to quasiparticle poisoning. 
Using Appendix B and mostly Eq. (|Blip . we value 



D. Small resistance 

We now turn our attention to the realistic situation of 
a small resistance so that cvi = Rk/Ri 3> 1. Of great 
interest to us is to understand how the quantum noise 
affects the long-time coherence of the EPR pair during 
the two-particle Breit-Wigner process involving the dots 
and the leads. An explicit calculation has been performed 
in Appendix D and we find an EPR current of the form 



I - = 0] 



exp(-27 e /ai) 



r(l + 2/ai 
where the crossed Andreev current at Ri 



I[Ri = 0] 



4e7g / sin(fc^^r) 



kp8r 



= readsi 

2<5rM. (32) 



\(p q \T>\DD)f « \T DL \* (er + e 2 - ir) f ^{e p + E cl /2) ^ 
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We have used the energy conservation e p + e q = — E c \ 
together with ei + e 2 = as well as e p + E cl /2 ~ 7;. 
Resorting to Appendix C, we obtain the following current 



I[5fi ~ E cl ] = 



4e 77l ( sm{k F 6r) \ 2Sr/< 



kpSr 



(30) 



Note that the Breit-Wigner resonance still occurs when 
ei + £2 ~ 2/ig ~ e p + e q reflecting the fact that 
for weak resistances the function Pi(E) diverges at 
E = and therefore the bath 1 only absorbs a tiny 
amount of energy during the EPR transportation process. 
The suppression factor (2/j, / E c i) 2Ri ^ Rk traduces the or- 
thogonality catastrophe arising in dissipative tunneling 
problems^iiiSiiiliSLill. Now, let us return to the discussion 
on the decoherence (decay) of a noisy EPR Cooper pair 
that is described, after averaging over the bath degrees 
of freedom, by the dot state (e -1 ^ 1 )^^^! 
In the weak resistance realm^i, we can evaluate 



Exploiting Eq. I|D8|) . we realize that compared to the 
noiseless case where R\ = R2 = the current becomes 



( e -i5<M*)) _ exp 



1 / £ el ait 
a.i \ ftfi' 



(33) 
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At very long-times t 3> hit / (pt\E c x) or very low energies 
(i <C E c i , by coupling to the environment the noisy EPR 
pair looses its phase coherence affecting drastically the 
efficiency of the two-particle Breit-Wigner resonance be- 
tween the dots and the leads LI and L2. It is worth to 
note the similitude with the power-law suppression in the 
Andreev entangler with Luttinger leadsi&ii. For larger 
resistances, at low fj,, we observe a precursor effect of the 
dynamical Coulomb gap mentioned above. 



It is important to visualize that in that case Pi 2 (E) yields 
a visible singularity at E = E c i (which is reminiscent of 
the situation where Ri is large and R 2 — 0) and therefore 
to get a current through the structure again one must re- 
adjust us — * Hs+E c i/2 and [i\ — > ni — E c i/2. We replace 
\(pq\T'\DD)\ 2 by its value in Eq. leading to 

/ = I[5n ~ E cl ] ( ' de q F 2 ( ei + e 2 - M; - e q ), (40) 



III. TWO INDEPENDENT BATHS 

We can straightforwardly generalize the previous anal- 
ysis to the case of two (independent) environments as- 
suming R 2 ^ 0. The two-particle Breit-Wigner transport 
between the leads and the dots gets modified as 

W DD j = | (pq\T'\DD)\ 2 Pi 2 (ei + e 2 - e p - e q ), (34) 



and therefore 



I[6n ~ E cl ,a 2 > 1] ~ I[6fj, ~ E cl ] 



where we have introduced 



dt 
2nh 



e iEt/n e Ji(t)+j 3 (t) 



(35) 



and J 2 {t) = ([6<fc(t) - Scj) 2 (0)}S<j) 2 (0)). When the two 
resistances are larger than Rk, P\ 2 {E) = 8(E—E c i—E c2 ) 
where E c2 — e 2 /(2C 2 ) ~ E c \ is the charging energy of 
the capacitor C 2 . To get a finite current between the 
superconductor and the leads 1 and 2 one needs to re- 
adjust 2fx s — > 2/i s + E cl + E c2 , //i — ► m - E c i, and 
fi 2 — > fi 2 — E c2 (we assume E c i ~ E c2 so that LI and L2 
are kept at the same electrochemical potential) . We find 



1 

E c2 



(36) 



We can observe a huge suppression factor ~ ( / y/E c i + 
l/E c2 ) 2 . When the two resistances are much smaller 
than Rk, which is the situation of most interest, we get 



PME) = 



exp(-27 e /a) 1 



T(2/a) 



E 



7T E 

a Ed 



2/ a 



where 



Ri_ 

Rk 



R2_ 

Rk' 



(37) 



(38) 



We obtain a result identical to that of a unique weakly- 
resistive bath with a.\ — > a (see Appendix D). Recall that 
the orthogonality catastrophe becomes more pronounced 
and the suppression factor in the crossed Andreev current 
now follows (2n/E c i) 2/a = {2fi/E c i) 2/oil+2/a2 ] consult 
Eq. (D15). Finally one could envision to investigate the 
asymmetric case where one resistance is well prominent, 
e.g., Ri, and the other is weak but nonzero resulting in 



exp(-27 e /q 2 ) 1 
r(2/a 2 ) E 



7T E - E c l 
a 2 Ed 



-| 2/a 2 



(39) 



exp(-27 e /a 2 ) 
r(l + 2/a 2 ) 



E-l ■ (41 » 



The crossed Andreev current becomes markedly sup- 
pressed by a factor proportional to (7 / E c i) 2 (2/j,/Eci) 2 ^ a2 
and ji <§; Ed is the applied bias voltage before the re- 
adjustment us —>■ i^s + Ed/2 and yu.; — ^ pt; — E c i/2. 



IV. PARASITIC DIRECT ANDREEV 
PROCESSES 

Thus far, we have completely omitted processes allow- 
ing the two electrons forming the Cooper pair to jump 
onto the same dot. In absence of thermal effects, the lat- 
ter are somehow reduced due to the traditional Coulomb 
blockade phenomenon on the dots as well as the super- 
conducting gap. Below, we will precisely discuss the 
quantum noise effect on those processes as well as the 
efficiency condition(s) of the noisy Andreev entangler. 

In the absence of quantum noise, direct Andreev pro- 
cesses (where the two electrons take the same dot) are 
suppressed by a factor {-fi/U) 2 and/or (ji/A) 2 compared 
to the crossed Andreev process as a consequence of the 
Coulomb blockade and the superconducting gap- Nev- 
ertheless, those direct Andreev processes do not suffer 
from a suppression resulting from the spatial separation 
of the quantum dots. We will use the terminology of 
Recher et ali by identifying two distinct direct Andreev 
processes: (I) In the first step, one electron tunnels from 
the superconductor to, say, dot 2, and in a second step 
the second electron also tunnels to dot 2. There are now 
two electrons on dot 2 which costs the Coulomb repulsion 
energy U ; this virtual state is suppressed by 1/U. Hence 
the two electrons leave dot 2 and tunnel to lead 2 (L2) 
one after the other. (II) There is a competing process 
that avoids double occupancy on dots but leaves an exci- 
tation on the superconductor that costs 1/A. Here, one 
electron tunnels to, say, dot 2 and then the same electron 
tunnels further into lead 2. Finally, the second electron 
tunnels from the superconductor via dot 2 into lead 2. 

We first concentrate on the tunneling process (II) in- 
cluding the quantum noise. The current from the 
superconductor to the final lead state takes the form 



Dp",f 



(42) 



p,p' 1 
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V 2 + SV 2 



FIG. 6: (Color online) Due to the Coulomb blockade direct 
Andreev parasitic processes are (diminished and) triggered by 
single-particle resonances between dot <-> lead. Noise will also 
affect the low-energy coherence of the single-particle current. 



We define the transition rate from \Dp") to |/) as 
2 5{e } + E s B -ei- e p „ - E% p ") and 



W {U) - \w (H > 

VV Dp"J ~ \ W Dp",f 



(II) 
' J Dp"J 



E 

71=0 



E (fB I (f\H' DL \Dp"a) \D BP ") {D B p" | (Dp"a\ (43) 



1 



ttI* 

- h dl; 



1 



irj -H ir)- H 



H' DL ) \Dp"a)\D B p" 



The final state |/) with two electrons in the same lead in 
the singlet state obeys |/) = (^/V2)(a p ^a p , ^ — a p ^a p ,^)\i) 
and we have introduced the intermediate state \Dp") = 
d)_ a a pll(J \i) with one electron on, say, dot 2, with spin — a 
and the other one already in lead 2 with spin a and mo- 
mentum p". The process where the first electron leaves 
the superconductor and tunnels to the lead 2 and where 
the second electron tunnels onto the dot 2 has to be 
accomplished in a very short-time ~ h/A and thus we 
assume that this process is instantaneous implying that 
the bath only reacts when the second electron resides on 
the dot 2. The intermediate state of the bath with one 
electron on the dot and the other in the lead has been de- 
noted \Dbp") and its energy Eg P . We also decompose 



(ii) 



i,Dp" 



(II) 
i,Dp" 



2 6(e p „ + e d +E 



Dp" 



Ci — E B ) where 



w 



(II) 

i,Dp" 



(D B p"\(Dp"a\ 



1 



ill - H 



■H'sd- 



1 



irj - H 



■H' DL (44) 



IT) 



Ho 



-H' 



SDI 



In the Lipmann-Schwinger operators, again we have 
taken into account that the bath's energies cancel out 



(in Hq now one must equate H\ ath — 0). From Sec. 
II B, we know that this is always justified in the weak- 
resistance limit which will be the situation of interest be- 
low. Note that a tunnel process from the state \i) to the 
state \Dp") does not have to be resummed further since 
this would lead either to a double occupancy of the dot 
that is suppressed by 1/U or to a state with two electrons 
simultaneously in the lead that is suppressed by a factor— 
7///i < 1. For convenience, we have suppressed the label 
I = 1,2 in the Hamiltonians H' D L and H' SD[ above. In 
the presence of voltage noise, we can yet decompose 



W 



(ii) 



Dp" J 



w 



(H)o 
Dp" J 



l \D BP ") 



(45) 



\ *U, \ E' B - <,-e p „-E°v") 



W {II)o P,(e „ 

vv Dp" ■,f Jr l\ t P" 



as well as 



(ii) 

i.Dp" 



W^ ) ; i \(Dgp"\e-^\ l g) 



(46) 



x <5(e P " +ei+ E B 



Dp" 



-ei- Eg) 
-ei), 



and (...)' // '° are exactly the transition rates occurring 
in the absence of dissipative effects. Note that the 
product of Hamiltonians H' SD H' DL H' SD is equivalent to 



-idipi 



HsdHdlHsd- Moreover we identify 



w (II)o (II)o 
W Dp"J W i,Dp" 



2^v s {T SD T DL y 
A(e p + ei - 
(e P + V)/2 



i-Yl/2) 
- £i - ill/2 



(47) 



+ £/ - «7i/ 2 



The momentum p" does not appear in the final result 
because one must satisfy either <5 PlP " or S PtP >. In the 
absence of noise where e p + e p > = 2fis = = e p " + 
ei we explicitly recover the formula of Ref. 0. Now, 
repeating the same calculation for the process (I) we find 
that WjyL, j — Wg*), f and furthermore w^g „ obeys 



w 



(I) 



%„ = (D B p"\(Dp"a\ 



-H'dl- 



irj -H ii] - Ho 

x-^H'sdIWb). 



H'sd (48) 



Compared to w^gpn the order of tunneling events has 

changed. W^g p „ has a form similar to Eq. H4bf> and the 

amplitude product Wgh tWfgZ,, is given by Eq. i|T7|) 
but with A being replaced by the Coulomb gap U/ir. 

In Appendix E, we have summarized the calculations 
of the direct Andreev current /' = I^ 1 ' + and the 
main result is that for relatively weak resistances 
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and for the noiseless case we again agree with Ref. 



i'[Ri = o] 



h£ 2 1 



(50) 



£~ l = 1/(ttA) + 1/U has been already mentioned in 
the introduction. It is important to stress that even 
though the direct Andreev current is suppressed by a 
factor (ji/£) 2 compared to the crossed Andreev current 
I, the former is less sensitive to voltage noise in the sense 
that this is only affected by an extra factor ~ (p/ E c i) 2 ^ a ' 
as opposed to ([j,/E c i) 2 / ai+2 / a2 for the crossed Andreev 
current / given explicitly in Eq. (D15) assuming two 
environments. This stems from the fact that the crossed 
Andreev current involves a two-particle Breit-Wigner res- 
onance and thus is more sensitive to voltage noise than 
the direct Andreev processes which demand that one elec- 
tron instantaneously leaves to the lead I hence producing 
a single-particle Breit-Wigner type transport through, 

e.g., Wpp), f. This is the uppermost issue of our paper. 

Note, in passing, that for the asymmetric case where 
i?2 = strictly whereas R\ would be finite (but much 
smaller than Rk), which means that dissipation only con- 
cerns electrons residing on dot 1, then we easily extract 



I'[ai > 1] 
I'[Ri = 0] 



cxp(-27 e /o: 1 ) / 2^i7r 



r(l + 2/e*i) \ ai E cl 



2/ai 



.(51) 



When the two electrons forming the injected Cooper pair 
take the dot 2, the direct Andreev current is equivalent 
to that of the noiseless case. Moreover the direct An- 
dreev current stemming from the passage of the two spin- 
entangled electrons through dot 1 is affected by the noise 
in the same manner as the crossed Andreev current. 



V. CONCLUSION 

In brief, in the setup of Fig. 1 we have thoroughly in- 
vestigated the effect of voltage noise produced by the elec- 



,25.33 



on 



trical circuits in the vicinity of the quantum dot 
the transportation of nonlocal charged-2e Cooper pairs 
(spin-based EPR pairs). We emphasize that even though 
electron spins in a semiconductor environment show un- 
usually long dephasing times approaching microseconds 
and can be transported phase-coherently over distances 
exceeding lOO/im. 30 , in the realistic dot-based Andreev 
entangler introduced in Ref. the voltage noise may af- 
fect the transportation of those EPR pairs through the 
charge degrees of freedom. Although the spin entangle- 
ment is preserved at long times, as a result of the entan- 
glement of the charge 2e with the electromagnetic noise 
the (noisy) Cooper pair object inevitably decays at long 
times. Assuming almost instantaneous tunneling events, 
we have been able to build a "P(E) theory" for this prob- 
lem along the lines of Ref. J27J More precisely, when 
investigating the Breit-Wigner resonances between the 
dots and the leads, we have kept only the dominant one 



"photon" contribution from each bath. For moderate and 
symmetric resistances R\ ~ R2 <C Rk, the condition for 
the noisy EPR entangler to be efficient [I /I' > 1) reads: 



Ri/Rk 

~ ) ( ) > fo6r). 



E 

E c i 



(52) 



Here, we have considered that 5r < 7r£ and two inde- 
pendent baths. It is important to recall that the crossed 
Andreev current is triggered by an "EPR-pair" Breit- 
Wigner resonance between the two leads and the two 
dots and is therefore (slightly) more sensitive to voltage 
noise than the parasitic direct Andreev processes which 
only involve a single-particle Breit-Wigner resonance be- 
tween, say, dot 1 and lead 1 (one electron has been in- 
stantaneously transmitted to the lead, e.g., due to the 
Coulomb blockade effect). Assuming that the Coulomb 
gap is large enough to satisfy this renormalized efficiency 
condition, it would be interesting to probe experimen- 
tally the long-time decoherence of the EPR pair through 
the orthogonality catastrophe factor (2/j,/E c i) 4R '/ Rk ap- 
pearing in the crossed Andreev current. We like to em- 
phasize that one could envision to exploit the dissipative 
GaAs heterostructures of Ref. |2?] to build a quantum 
dot in the proximity of a two-dimensional electron gas 
in low density which then serves as a tunable source of 
dissipation; interestingly, the resistance of the envion- 
ment can reach few fcfi. At this step, it is certainly 
important to also give our opinion on the case where 
the two dots would be subject to the voltage fluctua- 
tions of the same environment possessing a resistance 
R. We find that Eq. (19) would turn into Wddj — 

l( w |r|^)| 2 |(/ B | e 2 ^|c B )| 2 5( e/ - ei -6 2 + 4-^). 

It follows that the crossed Andreev current would be sub- 
ject to a more dramatic suppression ~ (2/i/ E c i) 8R / Rk 
whereas the direct Andreev current would exhibit the 
same power-law suppression as the two independent bath 
case. The efficiency condition of the entangler turns into 



Ed 



SR/Rk 



> (k F 5r) 2 



(53) 



It is relevant to note that the setup of interest to us 
is quite different from that of a BCS-superconductor di- 
rectly coupled to two highly-resistive normal leads be- 
ing described by two electromagnetic environments^. In 
that case, similar to Luttinger liquid leadsii*^, tunnel- 
ing of two spin-entangled electrons into the same lead 
is diminished compared to the crossed Andreev process 
where the pair splits and each electron tunnels into differ- 
ent leads. The reason is that when a charge 2e tunnels 
into the same lead I the tunneling process is accompa- 
nied by a phase e~ l2 ^ 1 ; this doubling of the phase leads 
to a more prominent suppression of current compared to 
the case where the two electrons take different leads. We 
insist on the fact that in the setup based on quantum 
dots, when the two electrons tunnel into the same lead 
the current cannot be triggered by charges 2e due to the 
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Coulomb blockade; the main process between the dot I 
and the lead I is a single-particle Breit-Wigner resonance. 

Note in passing that similar conclusions would typi- 
cally arise when including the noise in the normal leads 
LI and L2 in the setup of Fig. 1; more precisely, in that 
case we should carefully replace ai^ a — ► axka^ 5 ^ 1 ^ in Eq. 
(8) and hence the same conclusions could be derived. 

Since the spin entanglement is not really affected by 
the electrical noise, one could ask whether it would be 
possible to detect the spin entanglement of a noisy EPR 
pair despite the suppressed transmission probability of 
the latter at low voltage or long time. For example, Ref. 
IbH envisions to introduce a beam-splitter and focus on 
zero-frequency current correlations. Those quantities will 
be affected in a similar way as the currents (e.g., the total 
noise of the current flowing out of the superconductor 
is related to currents through the Schottky's result) 
and thus the detection of current correlations becomes 
highly dependent on the bias voltage. However, since the 
currents and the zero-frequency current correlations are 
affected in a similar way, we must admit that the Fano 
factors given in^] should not be modified that might give 
some hope to detect the nonlocal spin entanglement. 

Finally, the question whether the EPR pair can sur- 
vive if the dots are pushed away from resonance and are 
singly-occupied, is an interesting question that would be 
worthwhile to investigate further. The (related) struc- 
ture with a double dot in the Coulomb blockade regime 
coupled to two superconducting leads^i is well known to 
induce an antiferromagnetic coupling between the dots. 
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APPENDIX A: BCS NOTATIONS 

The s-wave superconductor is described by the BCS 
theory. The BCS Hamiltonian takes the form 

H s = Y, E *rL'r*» ( A1 ) 

kcr 



Ek = + A 2 being the quasiparticle energy and A 
the superconducting gap. The quasiparticle operator 7ko- 
is related to the electron annihilation and creation opera- 
tors Ckc and c k(T through the Bogoliubov transformation 

CkT = Mk7k| + «k7-kj,> ( A2 ) 
c-k| = Uk7-ki - Vk7kT> 

where the coherence factors itk = yl + (£k/-Ek)/v / 2 
and t>k = (£k/-E'k)/v / 2 have been introduced and 

£k = £k — MS is the normal state single-electron energy 



counted from the Fermi level /is. We choose energies 
such that fig = 0. Vv( r i) annihilates an electron in 
the superconductor at the site r/ and ij) a {r{) is related 
to Cka by the Fourier transform j/v(rj) = J2k elkTlc k<j- 
In our calculations, we will have to compute quanti- 
ties like (ilTko-V-aOri) = J2w e ' ik ' r ' (Al^aC^-a- The 
only terms which are non-zero should be proportional to 
(*|7ko-7kcr — Hence we infer that we shall select Ck' CT = 
c_kcr — > e ( ri'k7k CT m the equations (A2) above resulting in 
(i\-fkaCk'-a = (i|uk6(k + k')e CT (e a = ± for a =\,{). We 

evaluate (i\ip a (nhuc = Sk"(*l eik " r ' c k"<T7iL in a simnar 
way and we easily extract («|ck"fj7kcr — (*|wk<5(k — k"). 

APPENDIX B: CALCULATION OF (pq\T'\DD} 

Our aim is now to present a detailed calculation of 
(pq\T'\DD) without imposing e p + e q = = (equal- 
ity stemming from the energy conservation in the absence 
of dissipation). First, we can rigorously simplify 

( DD \Y,(: tT H *dl-. jtHdl) \DD) (Bl) 

1 

~ 1 - {DD\j^TT Hh Ll ^jr H DL \DD) ■ 

Hence, we exploit (DDI-. — -n- = (DD\-. — as well as 

(m(Hh lLl i^H DlLl +H* D2L27 ^H D2L2 )\DD) = 
\Tdl\ 2 Eim iri -^- ei resulting in 

(DD\^—H* DL ^—H DL \DD) = -. ,(B2) 

IT] —Hq IT) — Ho IT] — 61 — £2 

where we have introduced the self-energy 

E=\T DL \ 2 ^2(ir,-e l -e )c )- 1 . (B3) 

i,k 

Akin to Ref. 0, we can straightforwardly decompose E = 
5ReE — 17/2 where 7 = 71 + 72 and 3?eE ~ 7/ ln(e c //Lt;) 
can be neglected assuming that the renormalization of 
the energy level is small. This leads to the expression 

^\"7-#o in -H J 

1 _ ei + £2 - in 

~ 1 _ -i 7 /2 - ei + e 2 -il/2' 

irj — ei—e? 

Similar results hold for the one-particle resummation 

00 , ,n 

E • TT H h 2 L 2 h 7tHd 2 l 2 \pD) (B5) 

^ \tr)-H in -Ho J 

then providing 

00 / 1 1 \™ 

_ 1 _ e P + e 2 - ir\ 

- 1 _ -» 72 /2 - e +e2 _ ij2 / 2 - 
"J — Ep- £2 
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Again, \pD) stands for a\ pa d\_ a \i) and implies that the 
electron 1 is in lead 1 (LI) whereas the other electron 
resides on dot 2. We can also introduce the state \Dq) = 
^icr a 2q-crN) representing the state with one electron on 
dot 1 and the other one in lead 2 (L2) leading to 



,r H* DlLl -- TtHd^ \Dq) (B7) 

1 e q + ex - 17] 



1 __- ni /2_ e +£l -i 7l /2' 

WJ— 6 q — 61 



We can now proceed and compute 
(pq\T'\DD) = (pq\H DlLl \Dq) 



(B8) 



oo , 1 



-H 



1 



i7] -Ho DlLl ir)-Ho 



-H DlLl ) \Dq) 



x (Dg|- —H D2L2 \DD) 

17] - H 



rr "-DL • rr " 

i?7 -H IT] - Ho 



H DL \DD) 



{pq\H D2L2 \pD) 

oo , 

£ ( — ^ 



ir]-H " D2L2 irj -H 



Hd 2 l 2 ) \ P D) 



x (pZJU —H DlLl \DD) 

it] - Ho 



rr —DL~ rr - 

17) - H 17] - Ho 



H DL \DD). 



We obtain 

(pq\T'\DD) = T DL 



e q + ei - "7 



e q + ei - «7i/2 
7t>l \ ei + e 2 - «7 



(B9) 



?>/ - e x - e q / ex + £2 - *7/2 



e p + e 2 - «<7 / T DL \ e 1 +e 2 -i7] 



' ( P + £2 - «72/2 \ i7 7 - e 2 - e P / ex + e 2 - »7/2' 
and therefore 

<pg|T'|£>£>) = -T 2 DL ■ 1 , - 61+62 (BIO) 



' £ P + e 2 - 172/2 ex + e 2 - i"f/2 
1 ex + e 2 - 17] 



' e q + ex - «7i/ 2 ex + £2 - *7/2 ' 
This finally leads to 

(pq\T'\DD) = -Tl L 6 V f62 ". t? ;, x (Bll) 
ex + e 2 - «7/2 

e P + e q + ex + e 2 - i-f/2 

(e p + e 2 - i7 2 /2)(e q + ex - «7x/2) 

When neglecting the resistances of the leads that contain 
the capacitors C\ and C 2 , one can exploit that e q + e p = 



2]is = and recover the result from Ref. 

(pq\T'\DD) =-T 2 DL {e 1 +e 2 -i7]) x 
1 

(e p + e 2 - i7 2 /2)(e q + ex - Hi/ 2 )' 



(B12) 



APPENDIX C: CALCULATION OF 

([daidiT ± «Mu]T") 

Here, we would like to compute (l/v / 2)([d 2 xdx| ± 
d,2]dii\T n ) where the abbreviation (...) stands for 
This part is formally equivalent to 

{DD\T"\i) = -^(tUdaidit - da t d U ] x (CI) 
-H S d 1 . „ H S D 2 \i)- 



17] -Ho it]- Ho 

We can already evaluate {i\[d 2 id\-\ — d 2 ]d\i\ ^-Hq ~ 
(i\[d 2l d^ - d 2 ]d n ) iii _l i _ t2 leading to 

(DD\T"\i) = -L- x (C2) 



V2 it] - ex - e 2 



(i|[d 2 ;dxt - tfcTdiil-fffifDi- 



1 



17] - Ho 

Now, following Ref. 0, we insert a complete set of single- 
particle (virtual) states 



1 = ^lio d \- a \ i )( i \ d i-<y"fk<yi 



(C3) 



/k<r 



such that 
(DD|T"|i) 



1 



1 



^2 - ex - e 2 



(*][<feidi T -d 21 d n }H S D 1 (C4) 



^ZllA-a^^dl-al^a ^ H S D 2 \l) 



Now we can use (i\di-aJk^ ^pjj- = (^-^TW i7? _i k _ ei ; 
i?k being the energy of a BCS quasiparticle. Moreover, 
considering that the dots are at resonance we can ap- 
proximate i7] — Ek — ei ~ —Ek and obtain 

(DD\T"\i) = ~ (i\[d 2l d^ - d 21 d u ]H SDl (C5) 

V2 w? — ex — e 2 



Y rTE~^L d Lj i )( i \ d i-^^ H SD 2 \ 



Jker 



T, 



SD 



y/2 -IT] + ex + e 2 



i\[d 2i d x] - d 2] d n ] Y Y Y 



^ 4a' ( r l hkA-Ji) (*\dl-*irk*4a» ( r 2) l J > ■ 

The terms which survive are those with I — 2 and a' = 
—a" = a (—a is the spin polarization opposite to a): 

\{DD\T"\i)\ = ^ —Ih— (i\[d 2l d 17 -d n d n ] (C6) 
x Y ^ d L*-( r i)7L4- ff N>(*l7k^- CT (r 2 )d 2 - CT 4- ff N>- 
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Now we can develop <J/_ CT (r 2 ) = J2w e * k ' lc k'-er as 
well as \E , cr (ri) = ^ k „ e lk ri Ck",x, and exploit the 
precious equalities («|7kerCk'-£r = (i|vk<5(k + k')e CT and 
(i|ck"o-7k (T = (*|uk^(k — k") which have been demon- 
strated in Appendix A. Note that in agreement with the 
BCS theory we satisfy k' = — k". Hence we converge to 

\{DD\T"\i)\ = 4 , ?° , fll[<Mi t - d2 T d u ](C7) 
V2 + £1 + £2 



E^ k(ri - r2) ^U-.K>- 



kff 



Now we can resort to the important equality 

E<*l[<Mi T - d 21 d n }e a dl4_ a \i) = 2, (C8) 



then leading to 



\{DD\T h 



SD ^k^k p ik(n-r 2 ) 



-it] + ei + e 2 ^ -Ek 



(C9) 



In fact, another term which consists of exchanging the 
roles of Hsdi and Hsd 2 or r i an d r 2 in Eq. (CI) should 
be also included (see Eq. (12)). This means that 



\(DD\T"\i)\= ^ D Y (CIO) 



x Ukl ' k ^k(n-r 2 ) ( -,k,r, r. ; 



(Cll) 



and finally 



\(DD\T"\i)\-. 2x2T " 



SD ^™ cog(k(5r) (C12) 



17] + 61 + £2 4^ -Ek 



Now, we can exploit u k «k = ^g- V-^k — £k = A/(2E k ) 



|(fl£>lT"]i)| = 2 ^ T ^ ^^L C os(k.fr). (C13) 

' V 11/1 -MJ + £1 + £2 ^ 2El V ' V ' 

Here, we can assume that |k| rs fc^? and we note kp.<5r = 
kpSr sin 8 where Sr — \Sr\ such that 



\(DD\T 



»,..,, . 2 V2Tf^ L d0co S (k F 5rsin6). 



-17) + £l + £2 2A J 



(C14) 



We recover the crossed Andreev contribution of Ref. 

\(DD\T»\i>\K* , ^ Si ff\ (C15) 
V2 — + £i + £2 KFOr 

which is still valid in the presence of voltage noise. 



APPENDIX D: CURRENT CALCULATIONS 
FROM CROSSED ANDREEV REFLECTION 

1. No dissipation 

Without dissipation, from Appendices B and C we get 

(Dl) 



/ = ->'|T Di |^' Sin(W 



e \ -> 



kpSr 



1 



p + £2 - J72/2)(£q + 61 - «7l/2) 

Now, when £ p + £ q = we can rewrite 
1 

(e P + £2 - i72/2)(£ q + ei - ryi/2) 
1 f 1 



S{e p + £ q ). 



(D2) 



1 



61 + £2 - J7/2 \£ P + £2 - 172/2 
and hence the current turns into 

e 



-£ P + £1 - iji/2 



E 



2h< 



IT, 



DL\ 



77s 



(ei +£ 2 ) 2 + 7 2 /4 
1 



sm(kpdr) 



kp8r 



(D3) 



1 



p + £ 2 -i7 2 /2 -£ p + £1-171/2 

Owing to the Breit-Wigner resonance, we can eventually 
simplify £ p ~ e c ~ £1 and thus this results in 

9 



1 



1 



£ P + £2 - IJ2/2 
which results in: 



+ 61 - ryi/2 



= — 5(6 p - 6l),(D4) 

72 



77s 



sin(fci?5r) 



(D5) 



ft (£1 +£ 2 ) 2 +7 2 / 4 V k F^r 

This is the formula in Ref. 0. Exactly at the resonance 
condition for the dots where £1 = —£2 = the crossed 
Andreev reflection gives 



I = 



4e7g / sm(kpSr) 



h 7 



kpSr 



(D6) 



It is important to keep in mind that one requires the 
electrons residing simultaneously on the dots to be in the 
singlet state configuration so that results of Appendix 
C can be safely applied^. This is well satisfied when 
the distance between the dots < £. Indeed, the current 
carried by the S z — configuration of the triplet state 
on the dots would be zero because 



^(i|[d 2J di T + d 2 ^d ll ]e a d[ cr d : 



t I 

2 — CT I 



0. 



(D7) 



In the same spirit, to be non-zero Eq. (C7) assumes that 
the two injected electrons have anti-parallel spin configu- 
rations a" — —a' and opposite momenta k" = — k' which 
is well-satisfied for 5r < £. Therefore we can rewrite 



I[Ri = 0] = 



4e7| 



sm(kp5r) 
kpSr 



(D8) 
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2. Weak dissipation on dot 1 

For a weak resistance Pi so that a\ = Rk/Ri 3> 1 
the EPR current can be approximated as 



Pt(2/x f 



rfe q -Pi2(2M5 - ei - e q ), 

where the function Pi%{E) has been defined in Eq. (37). 
Assuming that P c x ~ E C 2 this leads to 



e q ) (p9)I[<*i »l]~Wi = 0] 



exp(— 2j e /a) ( 2\m 



T(l + 2/a) \aE cl 



2/a 



(D15) 



1 



Op + £2 - «72/2)(e q + ei - ryx/2) 



Since we consider that the dots are close to resonance, we 
have replaced P\{2^s — £ i ~ £ 2) ~ <5(2/^s — £x — 62) where 
the function Pi (P) (which yields a blatant singularity at 
E = 0) has been precisely defined in Eq. (|26|l . Now, we 
concentrate mainly on the effect of the quantum noise on 
the two particle Breit-Wigner resonance between the dots 
and the leads. For large a\ again the function Pi(P) has 
unambiguously a pronounced singularity at E = im- 
plying that e p + e q — > in the tunneling process between 
the dots and the leads. Using Eq. (D4) this allows us to 
approximate 



APPENDIX E: CURRENT DUE TO DIRECT 
ANDREEV PROCESSES 

Using Eq. 1471) . for weak resistances aj > 1 we find 
that the current I' = 1^ + /' Z U (j ue t direct Andreev 
processes (when the two electrons tunnel to the same dot) 
can be easily valued leading to 



2c 



EE- 

p,p' 1 

1 



1 2 s\Tdl\ 
£ 2 



■fl(2/i S -Cp-e P ')(El) 
1 



+ ei - in/2 



+ e; - iji/2 



2^|T DL | 4 7 I fsin(k F 6r) 



p,q 



(ex + £ 2 ) 2 +7 2 / 4 V k F$r 



Air 



x — <5(e p -ex)Px(2Ms 
72 



(D10) 



We have approximated P;(2/is — e P " —e;) ~ S(2fxs — £p" — 
e;) and we have explicitly introduced = l/(irA) + 
l/U. In the absence of noise, the energy conservation 
implies e p + e p > = ~ e; and therefore one can exploit 



We then converge to 



77s 



ft (£1 +£ 2 ) 2 +7 2 / 4 



sin(/cj75r) 
kp5r 



(Dll) 



de q Pi(2^ s - ex - e q ). 



We have used the fact that the dots are close to the res- 
onance condition ei + £2 = 0. Using Eq. (|26|) we obtain 



77s 



ft (£1 +£ 2 ) 2 +7 2 /4 V k F 8r 



sm(kp5r) 



(D12) 



exp(-27 e /ai) 



r(l + 2/ax) 
Hence we can summarize 



2/ai 



2/i 

Pel 



2/ai 



I[a x > 1] - /[Pi = 0] 



exp(-27 e /ai) / 2fnr 



r(l + 2/ai) Vai^cx 



2/ai 



(D13) 



3. Two independent environments 

In the case of two symmetric and moderate environ- 
ments such that Pi, P 2 <§; Rk, we get 



sin(A:j75r) 
ft (£1 +£ 2 ) 2 +7 2 /4 V 



77s 



(D14) 



£ P + £; - ill/ 2 
resulting in 



+ 



+ £; - ili/2 



4-7T 
ll 



£0, (E2) 



p,p' 



. 2 7 II^l| 
7i^ 2 



h us 2 



■6(e p + e p ,)8(e p - q) (E3) 
= 2e 7 |7 



Again, we reproduce the result announced in Ref. 0- In 
the limit of weak ohmic resistors Pi and P2 we estimate 



I'[Bi = 0] 

J [a; > 1] 2^ 



/'[P; = 0] 



E 



de q P(2//s-£/-£q) (E4) 

exp(-27 e /a;) / 2/i7r ^ 2/q ' 
r(l + 2/az) U/^c 



Note that for the most probable situation of symmet- 
ric environments, the suppression factor (2/i/P c /) 2 / Qi is 
less considerable than that for the EPR-pair current I. 
We argue that this is well justified in this setup because 
when two spin-entangled electrons tunnel onto the same 
dot the Coulomb blockade forbids charge-2e transport; 
this results in a single-particle Breit-Wigner resonance 
between, say, dot I, and lead I which is less affected by 
the baths than the two-particle Breit-Wigner resonance. 
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